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ABSTRACT 

We have obtained 1.64, 1.90 and 2.15 /zm narrow-band images of five T 
Tauri stars in the TW Hya Association (TWA) using the Near-Infrared Camera 
and Multiobject Spectrometer aboard the Hubble Space Telescope. Most of the 
T Tauri stars in our study show evidence of absorption by H 2 vapor in their 
atmospheres; in addition, the low-mass brown dwarf candidate, TWA 5B, is 
brighter at 1.9 fxm than predicted by cool star models that include the effects 
of H 2 vapor but neglect dust. We conclude that the effect of atmospheric 
dust on the opacity is important at 1.9 /xm for TWA 5B, the coolest object in 
our sample. The available evidence suggests that the TWA is 5-15 MY old. 
Comparison of the colors of TWA 5B with theoretical magnitudes as a function 
of age and mass then confirms previous claims that TWA 5B is substellar with 
a mass in the range 0.02-0.03 M . The accurate single-epoch astrometry of the 
relative positions and separation of TWA 5A and TWA 5B reported here should 
permit the direct measurement of the orbital motion of TWA 5B within only a 
few years. 

Subject headings: open clusters and associations: individual (TW Hydrae, 
CD— 33°7795) — stars: low-mass, brown dwarfs — stars: pre-main sequence 
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1. Introduction 

For more than a decade, the young star TW Hya has been an enigma since it lies in 
a region of sky apparently devoid of the raw materials to form stars, nearly 13° from the 
nearest dark cloud, yet it is unambiguously a classical T Tauri star QRucinski fc Krautter 



|1983|) surrounded by a great deal of cold dust ([Weintraub, Sandell fc Duncan 19891) an d gas 



( ^uckerman et al. 1995[ [Kastner et al. 1997 ). Recently, the TW Hya mystery was solved: 



TW Hya, along with other T Tauri stars found in an area of ~100 square degrees of the 
southern sky ( |de la Reza et al. 1989| , |Gregorio-Hetem et al. 1992| ), compose a uniquely close 



association of young stars known as the TW Hya Association ([Kastner et al. 1997| ). 

At a mean distance of only ~55 pc, the TW Hya Association (hereafter TWA) is almost 
three times closer than the next nearest known region of recent star formation. Given the 
likely age (~10 MY) of the TWA, these stars could harbor very young planetary systems 
with fully formed giant planets or low mass, brown dwarf companions, and may still be 
surrounding by circumstellar disks. In fact, there is substantial evidence for circumstellar 
gas and dust around several of these stars ([Weintraub, Sandell fc Duncan 1989| , [Zuckerman 



fc Becklin 1993| , |Zuckerman et al. 1995| , [Kastner et al. 1997|) . The relative proximity and the 



absence of significant interstellar or intra-molecular cloud extinction in the direction of the 
TWA make the prospects for detecting substellar companions around these nearest T Tauri 
stars much better than the prospects for similar searches for young low mass companions 
around T Tauri stars in Taurus- Auriga, Chamaeleon, Lupus or Ophiuchus, the next closest 
regions of star formation. In a recent study of the TWA, Webb et al. (1999) identified a 
total of at least 17 sources as members of the TWA. In addition, Lowrance et al. (1999) and 
Webb et al. have reported the discovery of a likely low mass brown dwarf companion (M 
~ 0.02 M Q ) to TWA 5A (= CD-33°7795) in a combination of ground-based and Hubble 
Space Telescope (HST) observations. TWA 5B is found almost 2" from TWA 5 A; thus, 
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despite its relative physical proximity (~100 AU) to the primary TWA 5B is amenable to 
spectroscopic and astrometric studies, uncontaminated by light from TWA 5A. 

In this paper, we report results from imaging the fields around five stars in the 
TWA, including the TWA 5 system, using the Near Infrared Camera and Multi-Object 
Spectrometer (NICMOS) and the HST. The goal of this program was to search for 
companions around these stars. Our choice of three narrowband filters centered at 1.64, 
1.90 and 2.15 /im was designed to enable us to identify cool and low surface gravity 
objects, including substellar mass companions, through their likely strong signatures of 
H 2 absorption at 1.9 /xm. 

2. Observations 

We obtained images of five star systems (Table 1) in the TWA using camera 1 (NIC1) 
and camera 2 (NIC2) of NICMOS between 1998 May 30 and July 12 (U.T.) Observations 
of each of the five stars were made identically. Using NIC1 and filter F164N, we imaged 
each target in a four position, spiral dither pattern, with an integration time per position of 
33.894 s. Three images were obtained at each position for a total integration time of 406.73 
s. We carried out identical observations using NIC1 and filter F190N, with an integration 
time per position of 43.864 s and a total integration time of 526.37 s. Switching to NIC2 and 
filter F215N, we again obtained three sets of four-position dithered image suites; however, 
for the F215N observations we changed the starting position for the dithered image suites 
in order to obtain a better median filtered image for subtraction of the thermal background. 
The integration time per position was 15.948 s for the NIC2 images, for a total integration 
time of 191.38 s. 
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3. Results 
3.1. Imaging 

We find no sources in any of our images other than the previously known five primaries 
and four secondaries, to limiting magnitudes of 18.3, 18.4 and 17.5 in the F164N, F190N 
and F215N images, respectively, at distances beyond ~1''3 at 1.64 and 1.90 fim and 2'' 3 
at 2.15 fim. The images of TWA 5 (Fig. 1) reveal how how easy it is to detect and image 
young, intermediate mass, brown dwarf companions around stars in the TWA, even without 
a coronagraph. In addition, all nine imaged objects appear as point sources (with FWHM 
of Q'.'U, 0"16, and 0"18 in the F164N, F190N and F215N images, respectively), with no 
evidence (after deconvolutions performed with point spread functions [PSFs] generated 
using the software package Tiny TimQ, direct subtractions of PSFs, and examinations of 
azimuthally averaged radial intensity profiles [Fig. 2]) of extended emission around any 
of them. Thus, although some of these stars appear to be surrounded by circumstellar 
material (e.g., TWA 1 is surrounded by a circumstellar disk of radius ~3" that is viewed 
nearly face-on; [Weinberger et al. 1999| , [Krist et al. 1999|) , we conclude that these direct, 



narrow band images are insufficiently sensitive to image circumstellar disks around these 
stars. 



4 http:/ /scivax. stsci.edu/ krist/tinytim.html 



- 6- 

3.2. Photometry 

We report our photometry for these observations in Table 1. The factors used to 
convert from NICMOS count rates to absolute fluxes and magnitudes^ were 5.376665 x 
1(T 5 Jy sec ADIT 1 for the F164N filter, 4.866353 x 1(T 5 Jy sec ADIT 1 for the F190N 
filter, and 3.974405 x 10~ 5 Jy sec ADU -1 for the F215N filter with zero point flux densities 
of 1033 Jy, 862 Jy, and 690 Jy, respectively. Photometry was obtained by measuring the 
total counts within a 0''5 radius aperture and then applying a correction factor of 1.15 to 
compensate for the flux that falls outside of this radiusQ. 

Figure 3 shows how the photometry of TWA stars through the J, H, and K broadband 
and the F164N, F190N and F215N narrowband filters relates to the near infrared spectral 
characteristics of late-type stars. Stars with lower effective temperatures have increasingly 
strong water absorption bands centered at 1.4 /im and 1.9 //m which are very effectively 
probed by this combination of broad and narrowband filters. 

We find no systematic differences between the K and F215N photometry (see Fig. 3, 
right panel), despite the slight difference in central wavelengths and large difference in 
bandwidth. On the other hand the stars are systematically brighter at F164N than at 
if -band, the most extreme case being TWA 5B. Finally, TWA 5B is the only star with an 
absolute flux that is clearly lower at 1.9 than at 2.15 /im. 

The left hand panel of Fig. 3 clearly shows how water absorption in the 1.35-1.55 and 
1.7 - 2.1 /im regions will strongly affect broad band H measurements but will have no effect 
on observations with the F164N filter (see also the library of near-infrared spectra published 

5 http:/ /www.stsci.edu/ftp/instrument_news/NICMOS/NICMOS_phot/keywords.html, 
version 1998, December 1 

6 |http: / / www.stsci.edu/ftp/instrument_news/JNKJMUS / nicmos_doc_phot.html 
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by Langon & Rocca-Volmerange 1992). Thus, the if -band and F164N observations reveal 
the presence of different amounts of water vapor absorption in the spectra of most of the 
stars in our sample. 

3.3. Astrometry 

For the four binary systems in our sample, we measured the intensity centroids for each 
binary component and transformed the cartesian positions on the array into offsets in Right 
Ascension and Declination of the secondaries from the primaries (Table 2). Except for TWA 
8, for which the binary separation is such that the companion only appeared in the NIC2 
images, the results presented in Table 2 are those obtained using only the NIC1 images 
since the spatial resolution is highest when using the NIC1 array. The NICl-based results 
in Table 2 are the statistical average and standard deviations based on measurements of the 
F164N and F190N images. The pixel to RA and Dec conversions were done using the plate 
scale measurement ephemeris generated by the NICMOS instrument teamQ. Comparison of 
the results for the F164N and F190N images indicate that, in most cases, we can determine 
image separations to an accuracy of 0.02 pixels (< 1 milli-arcsec). 

Because our astrometric results are obtained from unocculted HST images and with 
the highest resolution camera in NICMOS, these results are much more precise than offsets 
previously reported for these binaries. They are, however, consistent with previous results 
(Table 2). In the case of TWA 5B, Lowrance, Weinberger & Schneider (1999) recently 
independently determined that the offsets reported in Lowrance et al. (1999) and Webb et 
al. (1999) have a sign error in RA; the corrected values are reported in Table 2. 



http: / / www.stsci.edu/ftp/instrument_news/MlCMUS / nicmos_doc_platescale.html 



- 8- 



4. TWA 5B 

4.1. The age of the TWA 

We have constructed an H-R diagram for the TWA (Fig. 4) using the pre-main 
sequence tracks of Baraffe et al. (1998). This H-R diagram is quite similar to that presented 
by Webb et al. (1999), which is based on the pre-main sequence tracks of DAntona & 
Mazzitelli (1997); however, Fig. 4 appears to constrain the cluster age more tightly than 
does previous work on the TWA, presumably because of improved physics included in the 
Baraffe et al. tracks (see Baraffe et al. 1997 for a discussion). Specifically, virtually all of 
the stars, including TWA 5B, fall between the 3 and 10 MY isochrones. In comparison, the 
H-R diagrams of Webb et al. and Lowrance et al. (1999) indicate that the TWA stars have 
ages in the range 1-100 MY while Kastner et al. (1997) suggested that the likely age of the 
TWA stars is 10-30 MY, based on lithium studies (upper limit) and X-ray luminosities 
(lower limit). TWA 6 and TWA 9A, which lie together almost on the 30 MY isochrone, and 
TWA 9B, which falls near the 100 MY isochrone, are mild outliers in our and the Webb et 
al. HR diagrams and appear to be older than the other TWA stars (however, see Webb et 
al. for other possible explanations). 

What other information do we have to constrain the ages of the TWA stars? Soderblom 
et al. (1998) used the lithium abundance to place an age range of 5-20 MY and a most 
probable age of 10 ± 3 MY on TWA 4 (HD 98800; EW(Li A6708) = 0.36 A); Stauffer, 
Hartmann, & Barrado y Navascues (1995) used the strength of the Li line to assign an 
upper limit of 9-11 MY to TWA 11B (HR 4796B) while Jayawardhana et al. (1998) assigned 
an isochronal age of 8 ± 3 MY to this star; and Webb et al. (1999) measured similar Li 
EW strengths for 14 of the 17 stars identified as members of the TWA and, on this basis, 
suggested that they are all less than ~10 MY. The excellent agreement between the ages 
estimated from the Li EWs and those obtained from photometry and pre-main sequence 
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evolutionary tracks suggests that the age of the TWA is well constrained to be in the range 
5-15 MY. 

4.2. Mass and evolutionary status of TWA 5A and 5B 

TWA 5A is a M1.5 star (Webb et al. 1999) with T eff = 3700 ± 150 K (Leggett et al. 
1996). The distance to the TWA 5 system is presently unknown but can be estimated 
as 55 ± 9 pc from the measured parallaxes of four members of the association (Webb et 
al. 1999). The range of distances is consistent with the approximate angular dimension 
of the association. With d = 55 ± 9 pc and K = 6.8 ± 0.1, we find M K = 3.10 ± 0.41. 
By comparing these values of T eff and M K with the evolution sequences of Baraffe et al. 
(1998), we find M = 0.75 ± 0.15 M and an age of 2.5 to 6 MY for TWA 5A, assuming it 
is a single, pre- main sequence star (Fig. 4). On the other hand, Webb et al. report that 
TWA 5A is suspected to be a spectroscopic binary. If we assume that TWA 5A is binary 
with equal mass components, the mass of each component decreases to 0.7 ± 0.15 M and 
the age range becomes 6 to 18 MY. 

Since the Baraffe et al. (1998) sequence does not extend to substellar masses, we 
analyze the photometric measurements of TWA 5B with evolutionary models computed by 
Saumon & Burrows (unpublished). These models use the same interior physics as Saumon 
et al. (1996) and Burrows et al. (1997) with the distinction that the surface boundary 
condition is provided by the "NextGen" sequence of atmosphere models computed by 
Allard and Hauschildt for cool stars (Allard et al. 1996, Hauschildt, Allard & Baron 1999). 
The atmospheric structures provide a surface boundary condition for the interior models 
by giving a relation between the interior entropy (where the convective zone becomes 
essentially adiabatic at depth) and the surface parameters S(T e s,g). This relation plays a 
central role in controlling the evolution of fully convective stars. Colors are computed from 
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the synthetic spectra, and are therefore fully consistent with the evolution calculation. This 
evolution sequence was calculated for objects with solar compositions and masses between 
0.01 and 0.3 M , and is very similar to that of Baraffe et al. (1998) since it uses the same 
input physics (equation of state, atmosphere models, nuclear reaction screening factors, 
etc.). A limitation of the "NextGen" atmospheres is that they do not include dust opacity, 
which becomes significant for T cff ^ 2600 K. 

Figure 5 shows the evolution of the absolute magnitudes at /, J, H and K bands, from 
1 to 100 MY, based on the models of Saumon & Burrows (unpublished). Each curve shows 
the evolution for a fixed mass. The two dashed lines highlight the 0.02 and 0.03 M & models. 
The boxes show the photometric measurements for TWA 5B (Webb et al. 1999, Lowrance 
et al. 1999), with the height of the box representing the ±ler photometric error and the 
width showing the 5-15 MY estimated age of the association. The absolute magnitudes 
of TWA 5B assume a distance of 55 pc and the ±9 pc uncertainty is shown by the error 
bar in the upper right corner. All four bandpasses indicate that the mass of TWA 5B is 
between 0.02 and 0.03 M & with an upper limit of ~ 0.06 M Q if the TWA 5 system lies on 
the far side of the association and near the upper limit of our age estimate. Given this 
mass range and the estimated age of TWA 5B, the models indicate that its surface gravity 
is 3.8 < log#(cm/s 2 ) < 4.0. 

Stellar and substellar objects with M > 0.012 M Q undergo a phase of nearly constant 
luminosity which corresponds to the fusion of their primordial deuterium content (D'Antona 
& Mazzitelli 1985, Saumon et al. 1996). This phase lasts for 2-20 MY and contraction 

- with a consequent steady decrease in luminosity — resumes once the deuterium is 
exhausted. Figure 5 shows that TWA 5B is almost certainly in the deuterium burning 
phase of its evolution. 
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4.3. Colors of TWA 5B 

The UK colors of TWA 5B are consistent with its dM8.5-dM9 spectral classification 
based on 0.65-0.75 /xm spectra (Webb et al. 1999| , [Leggett, Allard fc Hauschildt 1998|) , 



and thus a temperature of T e g- = 2600 ± 150 K ( [Luhman, Liebert fc Rieke 1997] , [Leggett 



et al. 1996). In a J — H vs. H — K diagram, TWA 5B falls well outside of the observed 



sequence of very-low mass stars and brown dwarf candidates in the field ( [Leggett, Allard 



fc Hauschildt 199"5| ), while all other members of the association fall along the observed 



sequence of field stars. This indicates that the H magnitude for TWA 5B may be erroneous 
(by > ler) or that its relatively low surface gravity results in a redder H — K color. 

The narrowband infrared colors are shown in Fig. 6 along with the synthetic colors 
from the "NextGen" spectra. Each curve shows the colors for T e $ = 2600, 2800, 3000 and 
3200 K (from left to right) for a fixed gravity. The colors of TWA 5B are shown by the 
triangle with error bars. For the estimated T e g = 2600 ± 150 K and log g = 3.9 ± 0.1, there 
is a reasonable agreement for the F164N— F215N color but the models are ~ 0.4 magnitude 
too blue in F164N— F190N. Consequently, TWA 5B is brighter at 1.9 /zm than predicted by 
the models. 

The F190N bandpass falls in the middle of a strong H2O absorption band (Fig. 3) 
whose strength probably is overestimated by the "NextGen" models. Allard et al. (1997) 
compare a sequence of near-infrared spectra of late M dwarfs with their synthetic spectra. 
In all cases, the models overestimate the depth of the H 2 band, an effect which increases 
for later spectral types. While an inadequate H 2 opacity may be partly responsible 
for this effect, Tsuji, Ohnaka, & Aoki (1996) have shown that the condensation of dust 
in atmospheres of low T e s results in a source of continuum opacity which decreases the 
depth of the water absorption bands. New atmosphere models including dust opacity 
(Tsuji, Ohnaka, & Aoki; Leggett, Allard & Hauschildt 1998) indicate that its effects on 
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the spectrum (and on broadband colors) become discernible for T eff ^ 2800 K but remain 
moderate (~ 0.1 mag) at the effective temperature of TWA 5B (~ 2600 K). While current 
models including dust opacity may not fully account for the relatively high F190N flux of 
TWA 5B, the F164N— F190N color of TWA 5B is a strong indication of the presence of 
dust in its atmosphere. 

4.4. Astrometry of TWA 5B 

At a distance of 55 ± 9 pc, TWA 5B lies at a projected distance from TWA 5A of 108 
astronomical units. TWA 5A has a spectral type of Ml. 5 and a likely mass for the central 
binary of ~1.4 M while TWA 5B has an estimated spectral type of M8.5 ( [Webb et al. 
|1999|) and a likely mass of ~25 Mj up . Given this information about the TWA 5 system, the 
orbital period P of TWA 5B should be P ~ 1000 years. Thus, the angular motion of TWA 
5B, assuming a circular orbit of radius 1''96 viewed nearly pole-on, would be 0''013 yr _1 (or 
0''010 yr~ x if TWA 5A is a single star with a mass of 0.7 M and P = 1300 years). 

The NICMOS observations of TWA 5 were obtained on 25 April ( [Lowrance" et al. 19991) 
and 12 July 1998 (this paper), a difference of ~0.21 years. In only one-fifth of a year, the 
orbital motion of TWA 5B would have changed its position relative to TWA 5A by only 
~0''0027, too small for the positional difference to be measured at these two epochs. Thus, 
the differences between the positions we measured and the corrected positions from earlier 
epoch observations reported by Lowrance, Weinberger, & Schneider (1999) are strictly due 
to the relative accuracies of the different measurements. 

Although we have demonstrated that the positional change for TWA 5B is measurement 
error, not orbital motion, we also have shown that the position reported in this paper is 
a very accurate "starting" position for TWA 5B. In addition, our results show that it is 
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possible to measure the relative separation of these two objects to an accuracy of only a 
few thousandths of an arcsec. Thus, the orbital motion of TWA 5B should be measurable 
to a fairly high degree of accuracy with ground-based observing facilities equipped with 
adaptive optics or with the refurbished NICMOS camera. 

5. Summary 

To the sensitivity limits of these data, our images reveal no detectable circumstellar 
disks or infrared reflection nebulae, and no low mass stellar or substellar companions around 
stars in the five studied TWA systems other than the previously discovered TWA 5B. As 
for TWA 5B itself, our results suggest that this object has a mass in the range of 0.02-0.03 
M Q , in good agreement with the work of Lowrance et al. (1999). Finally, while our single 
epoch observations cannot demonstrate or measure the orbital motion of TWA 5B, they are 
more than accurate enough to permit the measurement of this motion, in combination with 
future epoch HST or adaptive optics, ground-based observations, with a baseline of only 
about a year. 

We thank the referee for thoughtful suggestions which improved the clarity of the 
manuscript, A. Burrows for computing the evolutionary sequences used in this work, and 
F. Allard, P.H. Hauschildt, I. Baraffe and G. Chabrier for making their synthetic spectra 
and models available. This research was supported by NSF grant AST93- 18970 and NASA 
grants NAG5-4988 and GO07861.01-96A and is based on observations obtained with the 
NASA/ESA Hubble Space Telescope at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS5-26555. 
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Table 1. Measured Photometry of TWA Stars 



TWA 


Common 


Spec. 


F164N 


H 


F190N 


F215N 


K 


Number 


Name 


TvDe a 


\ L1La BJ 


( ma? i° 




fmaf?^ 


( mat? l c 


1 


TW Hya 


K7 


7.50±0.01 


7.65±0.1 


7.53±0.01 


7.41±0.01 


7.37±0.07 


2A 


CD-29°8887A 


M0.5 


7.26±0.01 




7.39±0.01 


7.20±0.01 


7.18±0.07 


2B 


CD-29°8887B 


M2 


8.08±0.01 




8.12±0.01 


7.94±0.01 


7.99±0.07 


3A 


Hen 600A 


M3 


7.52±0.01 


7.60±0.1 


7.60±0.01 


7.37±0.01 


7.28±0.07 


3B 


Hen 600B 


M3.5 


7.90±0.01 


8.07±0.1 


7.98±0.01 


7.79±0.01 


7.80±0.07 


5A 


CD-33°7795A 


Ml. 5 


6.93±0.01 


7.06±0.1 


7.00±0.01 


6.86±0.01 


6.83±0.07 


5B 


CD-33°7795B 


M8.5 


11.76±0.01 


12.1±0.1 


11.83±0.01 


11.49±0.01 


11.5±0.07 


8A 


USNO 21A 


M2 


7.59±0.01 


7.72±0.1 


7.66±0.01 


7.53±0.01 


7.44±0.07 


8B 


USNO 21B 


M5 


d 


9.36±0.1 


d 


9.13±0.01 


9.01±0.07 



a from Webb et al. (1999). 
b this paper. 

c from Webb et al. (1999). Also, Lowrance, Weinberger, & Schneider 1999 report H = 7.2±0.1 for 

TWA 5A and H = 12.14±0.06 and K = 11.4±0.2 for TWA 5B. 

d TWA 8A and 8B are more widely spaced than the size of the NIC1 field of view; hence, at F164N 

and F190N we were able to obtain data only for the primary. 
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Table 2. Measured Binary Separations' 



Source ARA ADec Separation Position Ref. b 





(arcsec) 


(arcsec) 


(arcsec) 


Angle 




TWA 2B C 


+0.281+0.001 


+0.492+0.001 


0.567±0.001 


29.73°±0.03° 


1 




+0.3 ±0.1 


+0.5 ±0.1 


0.6 ±0.1 


31° ±3° 


2 


TWA 3B d 


-0.866+0.001 


-1.187±0.001 


1.469±0.001 


216.11°±0.03° 


1 




-0.8 ±0.1 


-1.2 ±0.1 


1.4 ±0.1 


214° ±3° 


2 


TWA 5B C 


-0.038+0.001 


+1.960+0.006 


1.960±0.006 


-1.11°±0.03° 


1 




-0.1 ±0.1 


+1.9 ±0.1 


1.9 ±0.1 


-3° ±3° 


3 




-0.04 ±0.01 


+ 1.95 ±0.01 


1.96 ±0.01 


-1.2° ±0.1° 


3 


TWA 8B f 


-1.220±0.014 


-13.162±0.019 


13.219±0.021 


185.30°±0.06° 


1 




-1.3 ±0.1 


-13.0 ±0.1 


13.0 ±0.1 


186° ±3° 


2 
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Table 2 — Continued 



Source 


ARA 

(arcsec) 


ADec 

(arcsec) 


Separation 
(arcsec) 


Position 
Angle 


Ref. b 





a from primary to secondary, based on centroid positions in NIC1 array images for TWA 2B, TWA 
3B and TWA 5B. As TWA 8B was only in the field of view in the NIC2 images, this offset is taken 
from those NIC2 images. 

b 1 = this paper. 2 = [Webb et al. 1999) . 3 = [Lowrance, Weinberger fc Schneider 1999| . Lowrance, 
Weinberger & Schneider report corrections of the position of TWA 5B originally reported in both 
Lowrance et al. (1999) and Webb et al. (1999). 

C NIC1 plate scale: X=0. 0431862 arcsec px _1 , Y=0. 0430120 arcsec px _1 , based on plate scale 
measurements from June 4, 1998. Observations obtained May 30, 1998. 

d NICl plate scale: X=0. 0431887 arcsec px _1 , Y=0. 0430144 arcsec px _1 , interpolated from plate 

scale measurements obtained on June 4 and August 6, 1998. Observations obtained July 1, 1998. 
°NIC1 plate scale: X=0. 0431897 arcsec px _1 , Y=0. 0430154 arcsec px _1 , interpolated from plate 

scale measurements obtained on June 4 and August 6, 1998. Observations obtained July 12, 1998. 

f NIC2 plate scale: X=0. 0759831 arcsec px _1 , Y=0. 0753005 arcsec px _1 , interpolated from plate 
scale measurements obtained on June 4 and August 6, 1998. Observations obtained July 9, 1998. 
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FIGURE CAPTIONS 

Fig. I. — NICMOS Images of TWA 5A/B, as seen through narrowband filters at 1.64 /im 
(F164N) (1 pixel = 0"043), 1.90 /im (F190N) (1 pixel = 0"043), and 2.15 fim (F215N) (1 
pixel = 0"076). All the small scale features seen in the images except for TWA 5B are seen 
identically in the model point spread functions for NICMOS. 

Fig. 2. — Azimuthally averaged radial intensity profiles for five of the TWA stars in our 
sample as seen in the a) F164N (1 pixel = 0"043), b) F190N (1 pixel = 0"043) and c) F215N 
(1 pixel = 0''076) filters. Each profile is normalized to the brightness of the central pixel in 
that profile. All profiles look identical to the PSFs generated by Tiny Tim, except for the 
presence of known companions. 

Fig. 3. — Photometry of TWA stars. Left panel: Synthetic spectra for a gravity of log g = 4 
and T cS = 4000 K (~K7) to 2600 K (~M6), from top to bottom, respectively (Allard et 
al. 1996, Hauschildt et al. 1999). Bars at the bottom show the bandpasses of the Johnson- 
Cousin J, H and K filters and of the narrowband NICMOS filters F164N, F190N and 
F215N. Right panel: Broadband (Webb et al. 1999, Lowrance et al. 1999) and narrowband 
photometry for a representative sample of TWA stars. The photometry is expressed in 
arbitrary flux units, normalized to the J flux. Boxes show the width of the broadband filters 
and ±1(7 error bars. Solid squares represent narrowband photometry results. The stars are 
ordered by spectral type, with later types at the bottom. 

Fig. 4.— Evolutionary tracks from Baraffe et al. (1998) and the stars of the TW Hya 
association. Solid lines show the evolution of stars of masses from 1.0 to 0.1 M Q in steps of 
0.1 M Q and for 0.08 M , from left to right, respectively. Isochrones for logt(yr) = 6.5, 7, 
7.5, and 8 are shown by dotted curves. The effective temperatures of the TWA stars were 
obtained from the Webb et al. (1999) spectral types and the Luhman & Rieke (1998) spectral 
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type-T efr relation. The error bars on M K reflect the uncertainties in the photometry (Webb 
et al. 1999) and on the distances. 

Fig. 5. — Evolution of M/, Mj, M H , and M K magnitudes for very low mass stars and 
brown dwarfs. From top to bottom, the curves correspond to masses of 0.2, 0.15, 0.125, 0.1, 
0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, and 0.01 M , respectively. Heavier curves show 
objects with M > 0.08 M Q which eventually become main sequence stars. 

Fig. 6. — NICMOS narrow-band color-color diagram. Colors calculated from the NextGen 
synthetic spectra (Hauschildt et al. 1999) are shown by the dots. Models with the same 
gravity are connected by a solid line, with logg(cgs) = 3.5 to 5.5 from top to bottom, 
respectively. The effective temperatures shown are 2600, 2800, 3000, and 3200 K, increasing 
from left to right. The observed colors of TWA5 B are shown by a triangle. 
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